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Abstract 

We discuss theoretical and phenomenological aspects of the use of bo¬ 
son propagators with energy-dependent widths in predictions for high-energy 
scattering processes. In general, gauge invariance is violated in such calcula¬ 
tions. We discuss several approaches to restore gauge invariance, necessary 
for a reliable result. The most promising method is the addition of the rele¬ 
vant parts of the fermionic corrections, which fulhlls all Ward identities. The 
numerical difference between this and other approaches is studied. A number 
of recommendations are given for LEP2 computations. 


1 Introduction 


Monte Carlo calculations for scattering processes for LEP2 and higher-energy 
colliders are required to have a precision of better than one percent. It is 
obvious that under such circumstances the assumption that the W and Z 
vector bosons are stable particles, produced on-shell, is a gross misrepresen¬ 
tation of the physics. Rather, one has to describe them as resonances, with 
a hnite width so as to avoid singularities inside the physical phase space. In 
held theory, such widths arise naturally from the imaginary parts of higher- 
order diagrams describing the boson self-energies, resummed to all orders. 
This procedure has been used with great success in the past: indeed, the 
resonance can be described to very high numerical accuracy. However, in 
doing a Dyson summation of self-energy graphs, we are singling out only a 
very limited subset of all the possible higher-order diagrams. It is therefore 
not surprising that one often ends up with a result that retains some gauge 
dependence. 

In itself, this is not necessarily a problem if one treads wearily. An exam¬ 
ple is the situation at LEPl. Here, a careful separation of the gauge-invariant 
subsets can be performed, leading to a result which has no signihcant gauge 
dependence. For processes that become important at LEP2, the situation is 
in several cases more complicated. Since gauge invariance is intimately con¬ 
nected with the high-energy behaviour of the theory, it is to be expected that 
the effects of gauge violation become worse if the scattering process under 
study contains a ratio of masses, or of momentum transfers, that becomes 
large. An example, which we shall study in this paper, is the production and 
hadronic decay of a single W in the process 

e’''e“ —*• e~UeW^ e~Ueud . (1) 

Here, the electron may emit a virtual photon, whose can be as small as m^, 
where rrie is the electron mass: with a total center of mass energy of ^/s avail¬ 
able, we have a mass ratio of s/m^ = O (10^^), large enough to amplify even 
a tiny gauge violation in a disastrous wayQ. An other, currently academic, 
situation, connected with SU(2) rather than U(l)e,m., is the gauge cancella¬ 
tion which prevents the cross-section for e’''e“ — W~^W~ from blowing up 
for high energies. 

In order to arrive at phenomenologically reliable predictions, various ap¬ 
proaches can be followed^. In the hrst place, we may try to convince ourselves 
that, for the particular problem under study, the situation is actually not so 

^This was noted already in Ref. ||i|, and investigated further in Ref. |^. 

^A few were investigated in Ref. 0 . 
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bad. For instance, this is the case in the above-mentioned LEPl processes. 
In processes like 

e+e" ^ 7 ,Z ^ /i+/r“ , (2) 

there is no obvions dangerons large ratio of masses at energies aronnd the 
Z mass, as the relevant ratio is s/M^. One might therefore hope that, by 
the imposition of a cnt on the electron scattering angle in the process (P, 
which effectively leads to a lower bonnd on the of the virtnal photon, the 
effects of gange violation can be mitigated. This is, for instance, an implicit 
assnmption made in the Excalibur Monte Carlo [Q. Of conrse, snch a hope 
has to be borne ont by comparison with a gange-invariant calculation. 

One may sidestep the problem by simply performing the calculation of 
the matrix elements without any width, and only at the end use some ad-hoc 
prescription like the following Let the mass and width of a boson be 

given by M and F, respectively, and its momentum by (F may depend on 

. Then, if we multiply the matrix element by — M^)/(g^ — M^ -|- iMT ), 
the pole at g^ = is softened into a resonance, at the expense of mistreating 

the non-resonant parts. It should be noted, that there are examples, where 
this ‘fudge-factor scheme’ leads to deviations up to 30% 0. 

Another way to sidestep the problem is to use the ‘hxed-width scheme’, 
i.e., to systematically replace l/(g^ — M^) by l/(g^ — -|- iMT), also for 

g^ < 0. This gives U(l)e.m.-current conservation, but it has no physical mo¬ 
tivation. In perturbation theory the propagator for space-like momenta does 
not develop an imaginary part. Moreover, the hxed-width approximation 
violates the SU(2)xU(l) Ward identities. Note, however, that this does not 
lead to a bad high-energy behaviour in 4 fermions, as the unitarity 

cancellations do not involve the masses of the W and Z bosons. In the case 
of —>• 6 fermions (e.g., Wl scattering) the occurrence of IF-mass de¬ 

pendent couplings means the unitarity cancellations are violated by a hxed 
width. 

A minimalist’s approach is to make use of the fact that the residue of the 
amplitude at the (complex) pole is gauge-invariant Bi- One can split the 
amplitude accordingly, and resum only this pole. In this way higher-order 
corrections can be included consistently |P]. However, this ‘pole scheme’ 
breaks down near thresholds, and has problems with the radiation of photons 
of energy ^ T. 

Finally, one may determine the minimal set of Feynman diagrams that 
is necessary to compensate for the gauge violation caused by the self-energy 
graphs, and try to include these. This is obviously the theoretically most 
satisfying solution, but it may cause an increase in the complexity of the ma¬ 
trix elements and a consequent slowing down of the numerical calculations. 
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For the vector bosons, the lowest-order widths are given by the imaginary 
parts of the fermion loops in the one-loop self-energies. It is therefore natnral 
to inclnde the other possible fermionic one-loop corrections |^. These 
fermionic contribntions form a gange-independent snbset and obey all Ward 
identities exactly, even with resnmmed propagators [12|. This implies that 


the high-energy and collinear limits are properly behaved. In contrast to all 
other schemes mentioned above, the scheme proposed here does not modify 
the theory by hand bnt selects an appropriate set of higher-order contribn¬ 
tions to restore gange invariance. 

To solve the problem of gange invariance related to the width, we only 
have to consider here the imaginary parts of these diagrams^. This scheme 
shonld work properly for all tree-level calcnlations involving resonant W- 
bosons and Z-bosons or other particles decaying exclnsively into fermions. 
For resonating particles decaying also into bosons, snch as the top qnark, 
gange independence is lost. For simplicity, we take all fermions in loops to 
be massless in the following. 

The jnstihcation, inclnding masses and the details of the proper resnm- 
mation and renormalization procednre, will be given in a later pnblication 
|T2| . The method has already been mentioned, and implemented in the s 
and t channel as a Monte-Carlo generator for the processes > ii'iqq', 

m 


in Ref. 


For the special case of qq' this approach was also nsed 


by Banr and Zeppenfeld 0, who fonnd that electromagnetic cnrrent con¬ 
servation was restored by fermion loops that essentially rescale the WW'y 
vertex. 

Althongh the proposed scheme is well-jnstihed in standard pertnrbation 
theory, it shonld be stressed that all reparation schemes are arbitrary to 
a greater or lesser extent: since the Dyson summation mnst necessarily be 
taken to all orders of pertnrbation theory, and we are not able to compnte 
the complete set of all Feynman diagrams to all orders, the varions schemes 
differ even if they lead to formally gange-invariant resnlts. It is then a nn- 
merical qnestion how mnch their predictions differ. 


The ontline of this paper is as follows. In the next section, we stndy 
the process of Eq. |l|, with emphasis on its small-angle behavionr. We show 
how gange invariance gets violated by the imposition of an energy-dependent 
width, leading to completely wrong resnlts. This is repaired by the inclnsion 
of fermionic corrections to the three-boson vertex. The electromagnetic cnr¬ 
rent is conserved again, and all Ward identities are satished. We discnss the 
connection between onr resnlt and that of Ref. 0. In section ^ we present 

^As the Ward identities are linear, we can separate the real and imaginary parts. 


3 


















numerical comparisons between the various reparation schemes for the pro¬ 
cess O- We hnish with a number of conclusions and recommendations. 


2 Gauge cancellations in e ^ e V(,ud 

In this section, we consider the process 

e"(pi) e+(fci) ^ e“(p 2 ) u{pu) d{pd) , (3) 


and especially concentrate on small scattering angles 9 for the electron. We 
keep the mass of the electron hnite, but shall neglect all other fermion masses 
(also that of the positron), so that we shall not have to worry about diagrams 
with Higgs ghosts connected to the positron or quark lines. The massive case 
can be treated analogously; this will be covered in Ref. |^. Under these 
assumptions, we have to consider the subset of four Feynman diagrams given 
in Fig. [^, which conserves the electromagnetic current. 





Figure 1: The four diagrams of the process e e"*" —»■ e Vf.ud which are con¬ 
sidered here. 

The matrix element M. is given by 

4 

M = Jf, , = , 

i=l 
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M^ = 4iQ^gl F^(pI) v_(ki)j^-^^^^-^j^v_(k2) U-(pu)70'^’-(Pd) , 

= -4iQ^g^ F^(p^_) u_(p^)j'^j^^^--^j^v_(pd) v_(ki)j0V.(k2) , 

= -4iQdgl F^{p^_) u.{j)u)i^ v_{ki)^pv_{k2) , 

= 4iglv_{ki)'^^v_{k2) u_{pu)'laV_{pd) ■ (4) 

The spinors are written in a compact form, U-{p) = 4(1 — 7 ^)m(p) and 

p+=Pu+Pd , P- = ki-k2 , q = Pi-p2 , 
[F^{s)]-^=s-Ml+i^^{s) , ( 5 ) 

where is the W mass and 7 ^ denotes the imaginary part of the inverse W 
propagator, which for the moment is introduced as a purely phenomenological 
device in order to avoid the singularities. The charged weak coupling constant 
g^ is given by gl = MlGp/V^, Qi is the electric charge of particle i, and 

= {pi-p2T^g''"^^+{p2-P3T"g^^^^ + {p3-piT^g^^^" ■ (6) 

If we use conservation of the charged current in the massless fermion lines, 
we may write 

W'’'‘(p+, -p_, -q) = (2p+ - qrg"-^ + - 2pj9'‘“ . (7) 

The photon source is given by 

JA* = 9lu{p2)Yu{pi) . ( 8 ) 

qz 

Note that the electrons can have two spin states each, but the massless 
fermions only contribute when they are left-handed. 

The matrix element, squared and averaged over the spins of the incoming 
fermions, reads 

(|A<7 = , 

^ [kp^ + pM + tH ■ 

spins Q. _ 

Note the occurrence of q~^: we may approximate 
, 2, mlR^ S-R^, 
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where '/S is the total energy, rrie the electron mass, 6 the electron scattering 


angle, and = {pu +Pd + ^ 2 )^: therefore, \q^\ can be as small as 0{ml). 


Let ns now consider the nnmerical behavionr of 
of Ref. 


Using the old approach 


14 , we dehne 


= Pi 


PP2 


/3 = pVp^2 


( 11 ) 


so that = 0 and (r)^ = (1 — therefore, as \q^\ becomes small, 

each individnal component of also becomes small. We may now write 




q: 


q\l-(5f 


2r'‘r>' + 


( 12 ) 

The nnwanted q ^ behavionr of the cross-section will be mitigated to a g ^ 
behavionr, provided 




(13) 


This conservation of electromagnetic cnrrent is seen to be extremely im¬ 
portant here: any small violation of it will be pnnished by a hnge factor 
O {S/ml). Mnltiplying q^ into the fonr diagrams of Eq. we obtain 


W = 

= Mo {(pI-pI)Q^ ^w(pI) ^w(P-) 

+ Qe pM+) - (Qd - Qu) pM-)} , 

Mo ^ M%g^^ . (14) 

By taking = 0, and considering the two poles at and at 

pl_ = , we get from the condition Eq. |T^ that 

Qw Qd Qu ) ( 1 ^) 


the obvions condition of charge conservation. Therefore, we have 

W = -i QeMo P^{pI) PwiP-) {'IwipD - 'JwiP-)) ■ (16) 

Cnrrent conservation is therefore violated nnless 7 ^y(p^) = 'y^{p‘i.). The 
most naive treatment of a Breit-Wigner resonance nses a fixed width approx¬ 
imation, with 

The nominal width of an on-shell W is given by 


T 


w 


E 

doublets 


Nf 


GfMI 

67rV2 


(18) 
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Figure 2: The extra fermionic diagrams needed to cancel the gauge-breaking 
terms. 

involving a sum over all massless fermion doublets with Nf (=1 or 3) colours. 
In this approximation, there is evidently no violation of electromagnetic cur¬ 
rent conservation. 

The difficulty with the hxed-width approximation is that it cannot be jus- 
tihed from held theory. Indeed, in held theory the propagator only develops 
a complex pole oh the real axis if we perform a Dyson summation of the self¬ 
energy graphs to all orders. This self-energy is inherently energy-dependent: 
to a good approximation^, we may write 

7w(s) = , s>0 , 

Iwi't) = 0 , t<0 . (19) 

Consequently, propagators with space-like momenta cannot acquire hnite 
widths in contradiction to the hxed-width scheme. 

The theoretically most satisfying way to restore gauge-invariance seems 
to be the addition of one-loop vertex-corrections, which cancel the imaginary 
part in the Ward identities. In the process above, this boils down to adding 
the imaginary parts of the diagrams of Fig. ^ These are given by 


M^ = 


IGtt 


p„{pl) PJpl) al Z iV; (Qj - <3„) Z”'’- , (20) 


where we included the appropriate colour factor for the doublet, Nf. Using 


"^Eq. W exactly takes into account the contributions of massless fermions, but it should 
be noted that above the W mass there is a contribution from the VEy self-energy diagram, 
which has to be treated perturbatively. 
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cutting rules, we calculate 


= dQTi 

ZiTT J 




(ri - g )2 


7 ^/ 27 " 


( 21 ) 


which is the imaginary part of the triangle insertions. The momenta ri and 
r 2 are the momenta of the cut fermion lines with = ri + r 2 . The expression 
satisfies the following three Ward identities 


2«/3(.p+ = 0 , 

Z“'’% = +j(pW-p 7 '‘“) ■ ( 22 ) 

Because of the anomaly cancellation we have no explicit contributions from 
the part containing 7 ®. Possible effects due to a top quark remain to be 
studied. Attaching the photon momentum to the sum of the diagrams 
Ad 5 gives 


hPadd 


2 

t Q.M, p„(pI) P„(pl) al Y. 

doublets ^ 

2 

iQM«PJpl)PAP-)r^-^ 
i QcMo P„(p\) P„(pt) (pI ) , 


(23) 


where we used the Ward identity of Eq. ^ the definition of the nominal 
W width, Eq. as well as the fact that the external charged currents 
are conserved for massless fermions. It is clear that the extra diagrams 
exactly cancel the imaginary part in Eq. Hence electromagnetic current 
conservation is restored. 

Using the fact that all external fermionic currents in process (^) are con¬ 
served one finds for the compensating correction (Eq. pl| ) 


Co 


+ g^^q'^ci -+- g^°‘q'^C2 + g'^^p+Cs , 


= 32 


Cl = 


pXp^^q^ 

A 2 




^^plp-q^ , 2 , 2,2 

10 -h p+ + p_ + g 


20p' 


A 

2 p- g 
A 


/o 


3 3 




/o + 4d_2l71 + l#j 

A ptq 3 ptq 


















C2 = -8 


C3 

/o 


= +8 




2 2 2 

p^ptq 


X 


2pl Pip> + 1 

A 

2q^ ^ + 1 


A 


f^ + 4P£P^ + ^P+-P- 
X 

/o + 4?!^ + 


3 piq^ q^ 

4 p+- q 1 1 


3 p^_q^ 


p 2 


2 2(p_- g) + yX' 


a/A \2{p- q) - x/Xj 


, A = 4(p_- qY — ^p'^^q^ 


( 24 ) 


This expression, inserted in Eq. ^ gives the correction to the WW'^ vertex 
to be used in explicit calculations. 

Now we discuss the result of Ref. |^. They computed the process gg' —>• 
^Vl'X| with dressed propagators for the IT’s and an on-shell photon. Note, that 
in the process of single IT production, one had g +p_ = p+. In the following 
process, we have p_ = g +p+. Hence, we have put g ^ —g with respect 
to the former definitions of the momenta. The invariant momentum squared 
flowing through both IT’s is positive and hence the running width is non-zero 
in both propagators. Without addition of extra diagrams, the corresponding 
amplitude will again not be gauge-invariant. Using the previous result in this 
section, it is easy to see that, with g^ = 0 , one gets for the two cut diagrams 
corresponding to the cut > 0 : 


3 a ^ ' 

\^p\p^ 


(. 


p+g -p+-qg'^ p+ 


16 P+P- 


3 [pig- - p_- g ) g" 


(25) 


with p+ > 0 , p^ > 0 , g^ = 0 and a = p\ — p^_. Note that the hrst term is 
proportional to the tree-level ITIT 7 vertex. The cut diagrams corresponding 
to the cut p^ > 0 are related by crossing symmetry. Adding the four cut 
diagrams, one ends up with; 


^ S 2 (g“^p(( + g^"g^ - g'^^g") . (26) 

cuts 3 


Inserting the overall factor and the fermion lines, one sees that the extra 
diagrams amount to a scaling of the ITIT 7 vertex with 1 + It 

should be noted, that the factorization of the correction is not universal. 

However, to get electromagnetic current conservation in the process (^, 
one can also effectively write the correction (p^) in this form. In the limit 
g^ —> 0, the overall factor multiplying the standard Yang-Mills vertex is then 
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given by 


(27) 


1 + i 


p\ — P- 


The parts in (|^) transverse to are dropped since they do not play a role in 
restoring electromagnetic cnrrent conservation. Only if one wonld allow for 
a negative running width in the f-channel, rather than taking 'j^{p‘^) = 0, 
multiplying the standard Yang-Mills vertex with an overall factor l+iT^ /M^ 
would give a result that respects electromagnetic gauge invariance. 

These simple, effectively factorizing prescriptions for restoring electro¬ 
magnetic gauge invariance may be easier to implement in a Monte Carlo 
generator. However, in general they violate the full SU(2)xU(l) gauge in¬ 
variance and, even more, upset the balance between the diagrams taking 
part in the unitarity cancellations at high energies. Hence, the validity is 
limited to the low energy range a/s = e.g., LEP2. In contrast, 

the factorized form obtained from Eq. |^, being exact, does not have this 
problem. 


3 Numerical results for e ^ e UgUd 

The process e“e+ —>■ e~UeUd has been studied numerically. The fermions are 
all taken to be massless, except for the electron, which has a mass me- The 

0.511 ■ 10"^ GeV , 

80.22 GeV , 

1/137.036 , 

1.16-10-® GeV^ , 

175 GeV , 

< no GeV . (28) 

The fermionic width of the W boson is computed using Eq. This gives 
T^ = 2.02773 ... GeV. 

The cross-section for e“e^ —*• e~ Veud for the different schemes for two 
values of the minimum electron scattering angle are given in table [^. 

Note that all schemes were computed using the same sample, so the dif¬ 
ferences are much more signihcant than the integration error suggests. One 
sees that in this case, once current conservation is restored the results for 
the total cross-section of the different methods agree to O From 

Fig. P it should be clear, that if we include running-width effects without 


input parameters are given below. 


me = 
= 


a(0) 

Gf 

50 GeV < 


10 



Scheme 

Gross-section [pb] 

f\ _ no r\ _ 1 no 

^min ^min 

Fixed width 

Running width, no correction 
Fudge factor, with running width 
Pole scheme, with running width 
Running width, with Eq. ^ 
Running width, with Eq. ^ 

.08887(8) .01660(3) 

60738(176) .01713(3) 

.08892(8) .01671(3) 

.08921(8) .01666(3) 

.08896(8) .01661(3) 

.08897(8) .01662(3) 


Table 1: Total cross-section for e e"*" —*• e Ugud in different schemes. 

taking into acconnt the correction of the Yang-Mills vertex, too many events 
are sampled for small valnes of q^. 


4 Conclusion 

The violations of gange invariance associated with a naive introdnction of 
a hnite width for nnstable particles can have disastrons conseqnences. We 
have indicated that, in the case of the vector bosons, this can be cnred in 
a fnlly consistent way by the inclnsion of appropriate fermionic corrections, 
e.g., to the three-vector-boson vertex. It has been shown explicitly in the 
case of massless fermions and the ITlTy vertex that the electromagnetic 
Ward identity is restored and cnrrent conservation holds. In the process 
—>■ e~9eud, in which gange-breaking terms are amplihed by 0(10^^), 
this is shown to lead to a correct resnlt. The differences between this scheme 
and other ways to obtain a gange-invariant result have been shown to be 
small, much less than T^/M^ in this specihc example. The correction to the 
WW^ vertex is given explicitly in Eq. for current conserving sources, and 
in a simplified factorized form suitable for this process at not too high energies 
in Eq. |^. These functions can be incorporated in other event generators for 
LEP2. 
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Figure 3: The effect of gauge-breaking terms in e e"*" —> e u^ud as a function 
of the virtuality of the photon. 




